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We study the thermal escape problem in the moderate-to-high and high damping regime of a system with a 
parabolic barrier. We present an equation which matches our numerical results accounting for finite barrier 
effects. We discuss some issues regarding the relation between escape rate and mean first passage time and 
apply our results to the study of force spectroscopy problems. 



I. INTRODUCTION 

Seventy years ago, Kramers^ proposed an equation 
for the thermal escape of a Brownian particle out of a 
metastable well. This problem is considered an issue of 
fundamental interest and it has received great attention 
from an impressive number of scientific areas^— . We 
present here a minor contribution to this vast field: a 
detailed numerical study, by using Langevin dynamics 
simulations, of the activation problem in the moderate- 
to-high and high damping regime of the system. 

This work is motivated by a recent numerical study 
of the Kramers problem at low damping^. There, im- 
portant finite barrier effects were described, a very slow 
convergence to the infinite barrier limit of the system was 
reported and the accuracy of the Drozdov-Hayashi the- 
ory£ was proved. As it is shown there, an accurate theory 
is of fundamental interest to understand experiments in 
this physical regime. 

On the other side of the same coin, the high damping 
limit of the problem is found. This is the typical case 
for many of the current friction^— and biological-physics 
worksi^— , where escape rate is the fundamental concept 
to understand force spectroscopy experiments. In many 
of these works the Kramers result for the escape rate at 
high damping: 




UgUb c -AU/k B T 
2lT"f 



(1.1) 



FIG. 1. Potential profile. Due to thermal fluctuations, a 
particle sited in state A is able to overcome a barrier AU and 
escape to state C. 



for the prcfactor and we recover the extensively used re- 
sult of Eq. (|1.1[) . The given results were obtained for a 
parabolic barrier— in the so-called infinite barrier regime 
of the system, where AU/ksT 3> 1. 

Another important concept in thermal activation the- 
ory is the mean first passage time (MFPT). This is de- 
fined as the mean time (t) for particles injected with zero 
velocity at point a, to reach a point b, where they are 
absorbed. In the ovcrdamped limit, given the potential 
U(x), the MFPT can be exactly computed as^ 



appears as the theoretical starting point to understand 
experimental results. 

In his famous worki, Kramers also derived an expres- 
sion for the escape rate in the moderate to high damping 
limit, 



(t) H n = 



dye u(y)/k B T / dze -u( z )/k B T^ 
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(1.3) 

There exist an important literature studying the re- 
lation between transition rates and mean first passage 
times^ ' 17 i 18 . It has been usually assumed, and theoret- 
ically proved, that there is a close relation between both 
quantities and in fact 



where the damping dependent prefactor is given by 

Ud(7M) = [1 + (7/2^) 2 ] 1/2 - (7/2w 6 ). There, co a 
and uJb are related with the potential curvature at the well 
and the barrier respectively, see Fig.[TJ AU stands for the 
potential barrier and 7 for the damping of the system. 
For high damping, j/cob 1, we have k KHD := uib/j 
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(1.4) 



whenever the injection and adsorbing points corre- 
spond to two potential minima separated by a barrier. 
Eqs. (|1.3[) and fj 1 .4[) allow for computing the exact es- 
cape rate in the overdamped limit of the system and for 



2 



comparing to the Kramers infinite barrier and other es- 
cape rate results. 

We devote section II of the paper to this end where a 
critical review of existing theories is presented. There, 
we propose a relatively simple equation, Eq. (|2.3[) . which 
allows for computing the escape rate for any barrier and 
damping in the moderate-to-high and high damping do- 
mains. Then, in section III we study the time to reach 
the barrier problem. This time, (t) m , differs from the 
previously calculated (t) by a factor that changes from 
1 to 2 when damping increases, and it is the relevant 
quantity in some situations. 

Armed with our previous results, we attack in section 
IV the so-called force spectroscopy problems. There, the 
system is continuously biased by an external field and the 
escape field probability distribution function is recorded. 
We compare results for the mean escape field based in 
different theories. We stress the differences between re- 
sults for field to reach the next metastable state and field 
to reach the barrier. We also show that damping times 
field ramp is the relevant parameter in the overdamped 
limit, and that there exists a value of this parameter be- 
yond which important noncquilibrium effects dominate 
the problem. We finish by showing results in the full 
damping spectrum. 



II. MODEL AND RESULTS 

We use Langevin dynamics simulations to compute the 
mean time for a particle with a given initial condition to 
reach a defined point beyond the barrier. To be precise, 
we will be first interested in computing the mean time 
for a particle sited with zero velocity in the minimum 
of a metastable well A to reach for first time the next 
available potential minimum C, see Fig. [1] For us this 
will be the definition of the relevant mean first passage 
time (MFPT) (t) of the problem. 

To be definite we study the dynamics of a Brownian 
particle in a metastable potential: 



mx + mjx 



av_ 

dx 



(2.1) 
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FIG. 2. (color online) Theoretical estimations of the escape 
rate in the high damping limit normalized by the Kramers 
result. Red line shows the exact mean first passage time result 
which accounts for the finite barrier effects in the system. We 
also show PT, Drozdov and exact TST results (see text and 
appendix). 



We compute the escape rate of a particle out of the 
metastable well (see Fig. |TJ) as a function of the damp- 
ing of the system and for different values of the 6 := 
AU /kgT ratio, which is the most relevant parameter of 
the problem. Some authors, using different approaches, 
have addressed the problem of including finite barrier 
effects in the theory 2 ^— . This is an interesting issue 
since in many cases escape occurs at relatively small 
barriers. In our model system, the energy barrier de- 
pends on the normalized external force / = F/Vq as 



AU = 2V 



1/2 

[l — / 2 ) — /cos -1 / . Another important 
parameter in the theory is the value of the curvature in 
the metastable minimum (w a ) and the barrier (w&). In 

this case uj a = tjf, = (Vo/m) 1 ^ 2 (l — f 2 ) 1 ^ ". 

We have performed molecular dynamics simulations of 
Langevin equation (12.111 . The equations were numer- 
ically integrated using the stochastic Runge-Kutta al- 
gorithm 2 ^ and the MFPT was computed typically over 
10000 realizations. As in our previous work^, we have 
used to = 0.35 and Vq = 0.155. In most of the simula- 
tions we use ksT — 0.01 and change F and 7^ 



where for the potential we use 



V(x) = V (l-cosx)-Fx, (2.2) 

and is the stochastic force describing the thermal 
fluctuations. Here we consider white thermal noise, 
(£(*)} = and (f (t)£(f )} = 2 mi k B T8{t - t'). 

The tilted sinusoidal potential describes a particle in a 
periodic potential subjected to an external field, a situa- 
tion of interest in many fields. In particular, Eq. (|2.ip is 
experimentally realized by a biased Josephson junction, 
and the high damping regime can be addressed with resis- 
tively shunted tunnel junctions or with superconducting- 
normal-superconducting onesi£. 



A. High damping 

First, we compare in Fig. [2] exact results from the 
MFPT theory [Eqs. (fLI)) and (fl~3"j) ] to the Kramers infi- 
nite barrier one, Eq. I|l.l[) . We show also the Pollak and 
Talkner— (PT) and DrozdoA* 2 ^ results, which include fi- 
nite barrier corrections to the escape rates [Eq. (|A5[) in 
the appendix]. As a reference we also plot the r^f* re- 
sult^ (see appendix). These results were proposed for 
the high damping limit of the system. The figure gives 
a quantitative estimation of the accuracy of the Kramers 
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FIG. 3. (color online) Numerical estimations of the escape 
rate in the moderate-to-high damping regime normalized by 
the Kramers result of Eq. (Q3J. Black line shows the exact 
mean first passage time result for the high damping case. 



high damping result for every barrier. As it can be seen, 
finite barrier effects are not very severe (below 10%) for 
normalized barriers above 5. However some care have to 
be taken when considering effects involving barriers be- 
low this value. Figure also shows the convergence of the 
escape rate to the Kramers infinite barrier result. 

As we also see in Fig. [3J PT theory gives a simple 
expression very close to the exact result except for small 
barriers. Drozdov gives an even better approach, however 
its final equation is more complex and only valid for very 
high values of the damping. PT approach fails at low 
barriers and Drozdov is a good approximation only in 
the overdamped limit. 
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B. Moderate-to-high damping 



FIG. 4. (color online) Escape rate at different 7 for 5 := 
AU/k B T =3, 5, 7 and 10, and prediction of Eq. (TOJ). (a) 
Escape rate as a function of 7. (b) Escape rate normalized 
by the high damping result as a function of 7/0;;,. (c) Escape 
rate normalized by the result of Eq. (|2.3p for AU/ksT = 3 
and 10 and prediction of the Kramers (dashed lines) and PT 
(dotted lines) theories. 



We propose below an expression for the escape rate 
which gives a good approximation for any barrier heigth 
and damping. Equation (|1.3[) is exact in the overdamped 
limit of the system. However there is no result for the 
MFPT at other values of the damping of the system. 

Our proposal is based in observing our numerical re- 
sults for the escape rate at different values of damping as 
a function of the normalized barrier AU/ksT. Figure [3] 
shows results for jr (inset) and r/r KMHD [note that, dif- 
ferent from Fig. [5J now we normalize the rate with re- 
spect to Eq. ([L2]) ] as a function of AU/k B T for T = 0.01 
and 7 = 0.1, 1, 10 and 100. There, we see first that our 
numerical results for high damping agree the theoretical 
prediction from the MFPT expression, Eq. (|1.3[) . 

Figure [3] also shows that, except for the 7 = 0.1 curve 
(we will show below that this value of damping corre- 
sponds to the moderate damping region where r(^) gets 
a maximum), all data approximately lie in a same curve 
when compared to the r KMHD prediction. Thus, we pro- 
pose the following simple extension of the MFPT results 



of Eq. (|1.3|) to the moderate-to-high damping region: 

7 \iHD = x k KMHD ( ) x — . (2.3) 

\u b J {t) HD 

The formula recovers the correct r HD — l/(t) HD result for 
7/W6 ^> 1 and incorporates the 7 dependence suggested 
in Fig. [3] As we will see, this expression gives account of 
our numerical results in the moderate-to-high damping 
regime of the problem in a simpler and better way than 
any previous approach. 

In order to evaluate the accuracy of Eq. (|2.3p , we have 
performed numerical simulations of the escape rate of the 
system as a function of 7 for different values of the bar- 
rier. In particular, Fig. [4] shows our numerical results 
obtained for four different values of the AU/ksT ratio 
(3, 5, 7 and 10). For better understanding, we plot re- 
sults in three different ways. Figure. Hla) shows the rate 
versus the damping. As we see, the theoretical expres- 
sion is valid for values of 7 larger than 0.1 approximately. 
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FIG. 5. (color online) (a) Theoretical (lines) and numerical 
(points) escape rates as a function of the normalized damping 
parameter for different values of the barrier over temperature 
ratio. Higher lines are theoretical predictions for AU/ksT = 
20. (b) Different theoretical approaches to the escape rate 
and comparison to our numerical results, for AU/ksT = 10. 
Red lines stands for Drozdov-Hayashi (low damping) and for 
our r MHD (high damping) results. Blue dotted lines stands 
for the BHL approach (low damping) and the KMHD one. 
Black line shows the simplest interpolation equation for the 
full damping range (see text). 



This value is at the lower limit of the moderate-to-high 
damping region. Smaller values of 7 sets in the moderate- 
to-low damping regime that we studied previously^. This 
kind of plots have been extensively used in the literature, 
however due to the log scale dependence of the rate it is 
unable to evaluate the precision of the theoretical equa- 
tion. 

Figure S^b) plots r/r HD as a function of j/ujb to check 
the validity of our ansatz in Eq. (|2.3D . Plotted in this way, 
all the curves approximately collapse in a single function 
given by (7/^) x k KMHD {^ /uj^). As expected deviations 
appears for the lowest values of damping. 

To finish, we present in Fig. |U[c) the direct comparison 
of our numerical results to Eq. (|2.3[) and PT and Droz- 
dov theories for 5 = 3 and S = 10. In this figure we plot 
the escape rate normalized by our theoretical expression, 



r/r MHD versus 7/W&. We conclude that Eq. (|2.3p repro- 
duces our numerical results in a better way than previous 
theories for any barrier and in the full modcratc-to-high 
damping regime. 

As we have seen, our numerical results are well fitted 
by Eq. (|2.3[) for damping values above 7 ~ 0.2 in our 
units. At this intermediate values of the damping pa- 
rameter we are into the so-called turnover region, which 
characterizes the transition between the moderate-to- 
low and the moderate-to-high damping regimes. There, 
the classical transition state theory establishes r TST = 
(w a /27r) exp-AU/k B T. 



C. Full damping 

In order to supplement the information given in previ- 
ous works, Fig.[5ja) shows the escape rate divided by the 
r TST result for a wide range of damping values with spe- 
cial emphasis in the turnover region. As expected, when 
plotted as a function of the normalized damping j/uib all 
curves collapse one each other only in the high damping 
limit. As barrier increases, the maximum moves to lower 
values of "f/tot- In any case rate approaches very slowly 
to the infinite barrier limit. The figure also shows the 
theoretical results studied in Ref. [fi and here. We see 
that the combination of the Drozdov-Hayashi^ result for 
moderate-to-low damping values and our extended mean 
first passage time result for the moderate-to-high ones, 
allows for a fairly good approximation to the numerically 
computed result for all damping of the system. 

Fig. (5Jb) shows the results for AU/ksT = 10 and 
also plots three other theoretical approaches. On the one 
hand we use the Kramers moderate to high damping re- 
sult and the Buttiker-Harris-Landaucr— (BHL) one. We 
see that this last deviates importantly in the moderate 
damping region of the system. To finish we also plot 
the result of the simplest analytical expression which ap- 
proach the escape rate for all damping, which is given by 
the following interpolation formulae: 

r -! = r - J 4. f- 1 (0 4) 

simple KLD KMHD ' \ ' J 

where for a linear-cubic barrier we have r KLD = 
(7.2j/2n){AU/k B T)c- AU / kBT and r KMHD is given by 
Eq. (|1.2[) . This is a simple expression which, as seen 
in the figure, gives a first glance good idea of the escape 
rate behavior in the full damping spectrum of the system. 



III. TIME TO REACH THE BARRIER 

The general expression for the MFPT depends on the 
injection, a, and absorption, 6, points. It is interesting 
to remark that the adequate choice for both numbers in 
order to keep Eq. (|1.4[) valid usually relies on relatively 
vague criteria since for particles starting in a metastable 
well far enough the barrier and adsorbed in another well 
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FIG. 6. (color online) Figure shows (t(a,b)) as a function of 
b for a = xa and & £ (a, a + 2ir). For comparison, results are 
normalized with respect to (t) = (t(xA, ic)) for each case. We 
set F — 0.05 and plot results for three different temperatures, 
resulting in different AU/ksT ratios, xa = sin _1 (_F'/Vb), 
barrier xb = tt — xa and xc = xa + 2n are marked by vertical 
segments. 



well beyond the barrier, the dependence on a and b is 
negligible. Thus a and 6 should be far from the barrier. 
Equation (|1.3[) allows for a direct study of the (t(a, b)) 
dependence. Figure [6] explores (t(a, b)) as a function of b 
for a = (the potential minimun position) and different 
values of the model parameters. As it can be seen, the 
MFPT changes importantly in the barrier region. As 
expected the time to reach the next potential minimum 
approximately doubles the time to get the maxima. A 
similar study can be done to know the dependence of the 
MFPT on the injection position of the system, a. 

There exist problems where the relevant physical mag- 
nitude is the mean time spent by the system to reach the 
maximum of the potential barrier (the transition state 
B in Fig. [T]). We call here this time (t) M to distinguish 
it from (t), defined above as the time to reach state C 
beyond the barrier. 

Let us introduce the factor e := (t) M / (t) for the ratio 
between both times. As can be seen in Fig. [51 a rough 
estimation in the high damping limit is that e HD ~ 0.5, 
showing that a particle in the maximum has a probability 
close to one half to cross the barrier and close to one half 
to move backwards to the original well. However, this is 
not exactly the case and as it is well knownii the so called 
stochastic scparatrix is sited beyond the barrier and dif- 
fers from the deterministic scparatrix. There, interesting 
dynamic effects appear in the Brownian dynamics of the 
syste m 27 ' 28 . 

The value of (t) M can be computed in the overdamped 
limit from Eq. (|1.3[) fixing absorption limit b at the poten- 
tial barrier position. We have numerically explored how 
e depends on the normalized barrier height AU /ksT and 
damping of the system. We should expect a weak depen- 
dence on the barrier at high enough values of AU/ksT. 
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FIG. 7. (color online) (t} M /(t) for different values of AU/k B T 
and 7. In (b) we also show (lines) the function proposed in 
Eq. {33]) for e HD = 0.5 and 0.545. 



With respect to the damping, e changes smoothly from a 
value close to 0.5 for high damping to 1 for small damp- 
ing. At small damping any particle that reaches the max- 
imum has enough energy to overcome the barrier and 
rapidly slides down to the next potential minima. 

Figure [TJa) plots e as a function of AU/ksT at differ- 
ent values of the damping 7. Figure [7|b) shows e as a 
function of the normalized damping 7/wb at different val- 
ues of the normalized barrier AU/kBT. We can see that 
e depends very weakly on the barrier except for values of 
AU < 3ksT. In which respects to the damping depen- 
dence, e changes from e HD ~ 0.55 at our parameters for 
7 > 10 to e ~ 1 for 7 < 0.1. 

In a similar way that for evaluating escape rates, it 
would be useful to get an analytical expression to com- 
puting the mean time to reach the maximum at different 
values of the parameters of the system. From Fig.JTJb) we 
see that for every barrier we can write (t) AI = / (7/075) (t). 
Following the ansatz used to derive Eq. (|2.3p we now pro- 
pose 
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FIG. 8. (color online) Mean force (F) as a function of the 
damping 7 for first passing the barrier top (red) and the bot- 
tom of the potential well (blue) in a force spectroscopy simu- 
lation at T = 0.01 and with force ramp 2.5xl0~ 7 . Lines are 
for theoretical predictions based in Eqs. (|2.3[) and (|3.1[1 . Top 
inset: comparison to different simple theories. Bottom inset: 
normalized barrier at the different mean escape force values. 



for the case of a parabolic barrier. There e HD = 
Whd/Wbd is computed from Eq. (|1.3[) . Figure [7] com- 
pares this equation to the numerical results for two pos- 
sible values of e HD , 0.5 and 0.545. Though agreement is 
worse than that for the escape rate case, to our knowl- 
edge this is the first attempt to propose an expression for 
the (t) M (7) dependence. 



IV. FORCE SPECTROSCOPY 

In a typical force spectroscopy experiment the prob- 
ability distribution function of the escape force P(F) is 
measured performing many experiments where force is 
continuously increased at a given rate F. From these 
results the mean escape force (F) and its standard devi- 
ation can be trivially computed. Such P(F) can be easily 
related to the escape rates r(F) as^fi 



P(F) = r(F) 



dF 
dt 



1 - 



P(u)du 



(4.1) 



Alike, escape rates can be computed from measured 
P{F). In this section we will show results of the mean 
escape force in our system for different parameter condi- 
tions and compare to available theories. 

As in the case of the escape rate and its relation with 
the MFPT, it is possible to introduce different definitions 
of escape force. In our simulations we start at zero force 
with particles in a metastable state (state A in Fig. [T]) 
and increase the force at a given ramp. Then we define 
the escape force as the mean force for which particles 
first reach the next metastable state (state C). However, 
in some circumstances the relevant escape force can be 



related with a different state, for instance the potential 
barrier (state B). In the first case we say that a particle 
has escaped if it has reached a new potential minimum, 
in the second if it has just passed the potential maximum 

Figure [8] shows the results of the mean escape force 
computed for both cases. As expected, (F) increases 
with the damping and both curves overlap for moder- 
ate damping and differ at larger ones. Figure also shows 
the theoretical results obtained from Eq. (|4.ip with r(F) 
computed as Eqs. (|2.3|) and (|3.1[) . in this last case using 
(-hd = 0.545. As expected, such curves give account of 
the numerical results in all the studied range. We want to 
remark that the difference between both curves at high 
damping is due to the factor close to 0.5 between the 
MFPT for reaching states B or C. Thus, a factor of 0.5 
in the prefactor of the rates have a measurable effect in 
the measured critical force value, see Fig. [5J This issue 
is not usually considered in the literature, as in many 
cases the prefactor is not carefully considered. In addi- 
tion, some force spectroscopy experiments are interested 
in fact in computing forces to reach the maximum in the 
potential landscape of the system, the so-called transition 
state, B. There, results are usually analyzed at the light 
of the Kramers high damping result, Eq. (|1.1|) . How- 
ever, according to our previous discussion, such expres- 
sion have to be corrected by a factor of 2. Such a factor 
is unwittingly integrated in uncertainties on the different 
parameters involved in the escape rate prefactor. 

We also compare our results with predictions from 
Eqs. (|1.2p and (jl.ip and analytical results by Garg2I, 
see appendix. Based on the small error found in the esti- 
mation of the system escape rates, we can advance that 
even the simplest theories give fairly good results. Such 
analysis is made in the top inset of Fig. [5] where the the- 
oretical results obtained from the Kramers high damping 
and the Kramers moderate-to-high damping results (red 
lines) and the first order Garg theory (black dashed lines) 
for high and moderate damping are shown. 

It is interesting to remark that in this numerical ex- 
periment escape happens at moderate to low values of 
the barrier over temperature ratio (see bottom inset in 
Fig. HI). Thus we see that results obtained from the infi- 
nite barrier limit give reasonable estimations for physical 
processes involving such small barrier values. 

Next, we are going to discuss an important issue, re- 
garding the validity of the available theories at very high 
damping values or force rates, where noncquilibrium ef- 
fects plays a role. At high damping escape events are 
very rare and, if the force is varied fast enough the crit- 
ical force of the model (for which barrier disappears) is 
reached before particles have reached the barrier. Then 
we can say that the escape problem is somehow ill- 
defined. 

These observations are studied in Fig. [S] where the high 
7_F limit is reached. At high damping the probability dis- 
tribution function P(F) depends on F and 7 through 
the 7F product. As seen in the figure, the nonequi- 
librium high jF regime appears in our simulations for 
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FIG. 9. (color online) Mean escape force as a function of "/F 
for different force ramps. Open and solid symbols state for 
two different escape force definition criteria (see text). Lines 
are theoretical predictions for both cases. 
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FIG. 10. (color online) Mean escape force as a function of 
7 for two different force ramps F = 1CP 8 (lower curve) and 
3.33 x 10 -7 (upper curve) and comparison with some theo- 
retical predictions. 



7 F > 1CT 4 ((F) > 0.13 and AU/k B T < 2). There we 
observe that an increasingly important fraction of the re- 
alizations do not escape before getting the critical force 
value (Vo = 0.155) for which the potential profile does 
not show a metastable well structure anymore. This 
situation is shown by the filled circles curves in Fig. 
where it is plotted the mean escape force of the fraction 
of particles which has escaped from the potential well for 
F < Vq = 0.155. An adequate theoretical calculation 
based on the MFPT results allows for reproducing the 
observed results even in this regime. To compare, we 
have also computed the mean force for the particles to 
advance from the starting point a to a point b = a + 2ir. 
Result is show by open circles in the figure and, as ex- 
pected, the theoretical results based in MFPT theory are 
also able to explain our results. 



For the sake of completeness we present a plot of the 
normalized mean escape force for a wide range of damp- 
ing values, 10 -5 < 7 < 10 3 and two different force rates. 
As shown before, the escape rate has a maximum for 
7 ~ 0.1 and thus the escape force presents a minimum for 
such damping value. For larger values of 7 our Eq. (|2.3[) , 
correctly estimates the escape force. For smaller values 
of 7, as shown ir£, the DH equation for low damping 
escape^ is valid. However, given the relative complexity 
of both approaches, we want to compare the results to 
those obtained with a simple interpolation formula, and 
the Garg analytical predictions. 

The simplest analytical expression which approach the 
escape rate for all damping is given by the interpolation 
formulae of Eq. (|2.4[) . As seen in the figure, such expres- 
sion allows for a very good approach to the computation 
of the mean escape force. In addition, figure also shows 
the simple approximate published by Garg (see also ap- 
pendix). We plot first order approach for moderate and 
high damping cases and second order for the low damp- 
ing one. Such results also give a good approach to our 
numerical results in a wide damping area. It is important 
to remark now that as ramp is decreased, the moderate 
damping area enlarges. Thus, the border between the 
different damping regions is also given by the force ramp 
and a good knowledge of both parameters is needed to 
correctly understand any experimental result. 



V. CONCLUSIONS AND DISCUSSION 

In this article we have presented a number of re- 
sults concerning the escape problem for high damping 
and its consequences in force spectroscopy experiments. 
Such problems include for instance cell to cell adhe- 
sion and molecular bond cxperiment a 13 ' 33 " — , dissociation 
of molecular complexes 3 ^, mechanical unfolding of pro- 
teins^, single-molecule pulling experiment o 14 ' 37 " — , free- 
energy reconstruction 4 ^—, and study of friction at the 
nanoscaleSrildZ-^L^. 

Our work gives an accurate expression for escape rate 
in the moderate-to-high and high damping domains and 
allows for a quantitative evaluation of errors made in the 
use of different approaches to the problem. Infinite bar- 
rier results are usually assumed in literature without an 
estimation of the real normalized barrier values involved 
in any case. We have explored this issue and showed 
that finite barrier effects do not importantly affects force 
spectroscopy results though escape indeed may occur at 
very low barrier values (Fig. [8]) . This was also obtained 
in Ref. M- 

Another important issue is the role played by the 
damping parameter. This is a pretty difficult issue to 
evaluate in many real systems, although the validity of 
high damping approach is usually considered. We see 
that, on the one hand, high damping theories works 
fine for dimensionless damping j/u> > 2. On the other 
hand, if the jF product goes above 10 -4 , important non- 
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equilibrium effects appear and usual escape rate theories 
are not valid any more. We have pointed out the exis- 
tence of a high 7_F limit beyond escape problem is not 
well defined. 

Frequently, the main source of uncertainty in the prob- 
lem is the potential profile itself. Then, sometimes it is 
assumed r ~ roexp(— AU/ksT) with tq some empirical 
prefactor. In such cases, all interest is focused on the bar- 
rier value and the field dependence of this barrier, and 
the weak r$ dependence with F is neglected. However, in 
general, Vq contains information about the damping, the 
bias field and the barrier shape. 

In adhesion and other similar problems, free energy 
is frequently modeled by a tilted Morse-type potential. 
Then, the applied field creates a barrier—. However, 
next potential minimun is sited at +00 and escape rate 
is vaguely defined in terms of being far (but not too far) 
from the potential barrier. Then, mean time to reach the 
maximun can be an useful concept and sometimes, when 
a cusplike barrier is assumend, it is the only meaningful 

14 

We have discussed the relation between escape rate and 
MFPT and computed results for MFPT to the barrier. 
We argue that in some cases such parameter should be 
the relevant one. In close connection to this, in force 
friction experiments for instance, mean friction force is 
computed from the maximun force in each stic k 8 ' 9 ' 51 ' 52 
slip cycle. To compute such a force, the MFPT to reach 
the maximun is the relevant physical magnitude. This 
has been overlooked in the past, introducing an error by 
a factor of 2 in the calculus which can be easily added to 
the uncertainty in the evaluation of the system effective 
damping. This has been also the case in some analysis of 
single-molecule force experiments 3 -^. 

Most of the previously discussed issues plays also an 
important role when analyzing single-molecule pulling 
experiment a 37 ' 40 ' 54 . An interesting question here is that 
of the reconstruction of free-energy landscapes 4 ^. More- 
over, in such, and other systems, additional issues beyond 
the scope of this manuscript as diffusion in a multidimen- 
sional landscap o 41 ' 53 plays also an important role. 



1. Pollak and Talkner. 

Pollak and Talkner proposed 2 ^ the following expression 
for the escape rate: 



, fix) 

' PT — K KMHD X 1 TST 1 1 



(Al) 



where x = AlI/hsT and x t ne following function of 7 

7M 



(1+ ^2 /4) l/2 
7/2 



(A2) 



and 



fix) 



1 

3% 



2-3 X -i(x + l) 3 + 



3 (X-1) 2 (X + 1)(3X 2 + 12X+1) 
2 9x 2 - 1 



(A3) 



In the high damping limit \ ~> 1 an d fix) ~* —5/36 
then 



(A4) 



In this theory it is necessary to compute the term 



dx e-W*) 



(A5) 

Note that in the high barrier limit r™ 3 ^' ~ 
77*- c~ AU / kBT := r TST , the usual result given by the tran- 
sition state theory and used in Eq. (|1.2j) . 



2. Drozdov 



The problem was also studied by Drozdov in a series 
of papers^ ' 22 ' 32 where he proposed: 



r D — K D X r TST ■ 



(A6) 
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being n = 8/7 a good choice, 



L0 e = V^Tr//? 



+00 



-1/4 



dx e 0v W 



(A7) 



(A8) 



Appendix A 



and 



Here we will give some numerical expressions that have 
been used along the manuscript. 



, ,2 o3/2 f+00 

1 dx V' 2 {x) e^^\ (A9) 



2tt 



9 



0.14 
0.12 
0.1 
0.08 
0.06 



10'" 10' 3 10' 2 10~ 1 10° 10 1 10 2 

7 

0.06 





0.00575 
0.0055 

0.00525 
0.005 

0.00475 
0.0045 





.... . 









10"' 10' 3 10" 2 10' 1 10° 10' 10' 
Y 



FIG. 11. (color online) Evaluation of Garg results for two 
different force ramps. Points stand by the exact numerical 
evaluation, solid lines for first order approximation, dashed 
lines for second order approximation results. 



3. Garg 

To finish, we will present a short comparative study 
of Garg^i predictions for the mean escape force and 
its variance. There, both quantities are given as a 
series expansion of a certain parameter. Figure 1111 
shows a numerical study of the results reported by 
Garg for two different force ramps. Barrier height is 
given by 4v / 2Vb/3(l — F/Vq) 3 ^ 2 , oscillation frequency by 

2 1 / 4 (VVm) 1 / 2 (l-F/y ) 1/4 , T = °- 01 and ^ = 2-5 x 10- 7 
and 10~ 12 . Figure shows that for the case of (F) first 
order approximation is very good for the moderate and 
high damping cases (even better that considering the sec- 
ond order approximation). On the contrary, in the low 
damping result only the second order approximation cor- 
rectly estimates the exact escape rate. With respect to 
the variance, result are not that good. For the case of 
moderate and low damping first order approximation is 
not that bad and second order is better. However, for low 
damping first order result is not good but second order 
is much worse. 
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